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Abstract 

We consider Bondi's radiating metric in the context of the telepar- 
allel equivalent of general relativity (TEGR). This metric describes the 
asymptotic form of a radiating solution of Einstein's equations. The 
total gravitational energy for this solution can be calculated by means 
of pseudo-tensors in the static case. In the nonstatic case, Bondi de- 
fines the mass aspect m(u), which describes the mass of an isolated 
system. In this paper we express Bondi's solution in asymptotically 
spherical 3+1 coordinates, not in radiation coordinates, and obtain 
Bondi's energy in the static limit by means of the expression for the 
gravitational energy in the framework of the TEGR. We can either 
obtain the total energy or the energy inside a large (but finite) por- 
tion of a three-dimensional spacelike hypersurface, whose boundary is 
far from the source. 
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I. Introduction 

The concept of energy in general relativity is considerably more intricate 
than in any other branch of physics. Any physical phenomena, except grav- 
itation, is defined and described on a specific space-time, which is usually 
the flat space-time. For these phenomena the concept of energy can be in- 
tuitively conceived and mathematically realized. Generically, energy is an 
attribute of some physical system whose dynamics takes place on the space- 
time. Gravitation, however, acquires a distinct status because the dynamics 
of the gravitational field is the dynamics of the space-time itself. Conse- 
quently, the definition of the gravitational energy is not straightforward. 

The several attempts at defining the gravitational energy (pseudotensors, 
quasi-local energy, actions and Hamiltonians with surface terms) all agree in 
predicting the total energy of asymptotically flat gravitational fields. More- 
over, there seems to exist a predominant point of view according to which the 
gravitational energy is not localizable, i.e., that there does not exist a grav- 
itational energy density. These are probably the only two features shared 
by the various approaches, which are mostly based on the metric tensor. 
However, the very concept of a black hole lends support to the idea that 
gravitational energy is localizable. There is no process by means of which 
the gravitational mass inside a black hole can be made to vanish. 

A detailed analysis of the structure of the pseudotensors shows that the 
(covariant) gravitational energy-momentum tensor would have to be defined 
by means of the first derivative of the metric tensor. But it is well known that 
it is not possible to write down a non-trivial covariant expression involving 
the first derivative of the metric, which captures the energy content of the 
field. However, it is possible to write down such covariant expressions with 
tetrads, and M0ller noticed this fact long ago|], H, ||. 

The question of localizability of the gravitational energy can be dis- 
cussed in the framework of the teleparallel equivalent of general relativity 
(TEGR)0, ^, |5], 0], which is an alternative geometrical formulation of Ein- 
stein's general relativity. The action integral of the TEGR is constructed 
entirely out of the torsion tensor. The analysis of the canonical structure of 
the TEGR0 indicates the existence of a perfectly well defined gravitational 
energy density. Such existence is possible in principle, because the TEGR is 
defined in terms of tetrad fields. The torsion tensor allows the construction 
of a total divergence that transforms as a scalar (energy) density. In the 3+1 
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formulation of the TEGR the integral form of the Hamiltonian constraint 
equation C = can be written as an energy equation of the type|| 



C = H - E = , 
where E is the gravitational energy defined by 

where e = det(e^)i), {e(k)i} are triads restricted to a three-dimensional space- 
like hypersurface S, and T % is the trace of the torsion tensor: T % = g lk T k = 
gik e (i)j T(iy k , T(i)j k = dje^k — c^e^j. V is an arbitrary three-dimensional 
volume of integration and G is the gravitational constant. This expression 
is simple and powerful. It has been successfully applied to rotating black 
holes ||, de Sitter space (JToJ and conical space-times Ullj . The definition of 



gravitational energy in the TEGR may not be intuitively clear, but it is sup- 
ported by its mathematical simplicity and by the applications to the space- 
times listed above. The use of (1) requires only the construction of the triads 
{ e (k)i} with the appropriate boundary conditions, and which transform under 
the global SO (3) group. 

The torsion tensor that appears in the Hamiltonian formulation of the 
TEGR is related to the antisymmetric component of the connection T l - k = 
e^ m i % djet m )ki whose curvature tensor is identically vanishing. Such connection 
defines a space with teleparalellism, or absolute paralelism, or else fernpar- 
allelismus, according to Schouten||T2| . 

In this paper we investigate the energy of asymptotically flat gravitational 
waves, described by Bondi's radiating metric [13|]. Since the metric describes 



an isolated system the application of (1) is possible as it stands, provided we 
consider the metric in the 3+1 spherical coordinates (t, r, 9, 0) at spacelike 
infinity, for which t = constant defines a spacelike hypersurface. We note 
that the use of cartesian (rectangular) coordinates in the asymptotic limit is 
necessary for the evaluation of pseudotensors out of this metric. 

We recall that the Arnowitt-Deser-Misner (ADM) energy p^j is not suit- 



able for the analysis of gravitational radiation because it gives the total en- 
ergy of the space-time, both from the source and from the emitted radiation, 
whereas the Bondi energy evaluated at null infinity furnishes only the energy 



2 



of the source, from which it is possible to derive the well known formula for 
the loss of mass. 

The relevance of the definition (1) resides precisely in the fact that we can 
evaluate it on a large but finite volume V of the three-dimensional spacelike 
hypersurface, thereby not including the emitted radiation outside V. In view 
of the field equations (which are not considered here), the energy inside V 
turns out to be a decreasing function of time. 

It is important to remark at this point that Bondi energy has been cal- 



culated in several geometrical frameworks, by different aproaches|15|, [16], [17 



|T8| , [HJ. A commom feature of these approaches is that they yield the total 
energy of the field. In contrast, we will consider finite volumes of spacelike 
surfaces and obtain the energy contained within large spherical surfaces of 
radius r Q up to the — term. 

In the next section we briefly describe the Lagrangian and Hamiltonian 
formulations of the TEGR. In section III we compare our energy expression 
with M0ller's expression. We show that both expressions agree for the total 
gravitational energy, but inspite of similarities they disagree when applied 
to finite volumes of the three-dimensional space. In section IV we write 
Bondi's metric in (t, r, 9, 0) coordinates at infinity and proceed to carry out 
the construction of triads for the spacelike hypersurfaces S. There exists 
an infinit number of triads that lead to the metric restricted to the three- 
dimensional hypersurface. However, only two of them will be considered 
in detail. In section V we calculate both the total energy of the field and 
the energy contained within a large sphere of radius r . The total energy 
obtained by means of (1), in which case the integration is made over the 
whole S, agrees with the known result for the Bondi energy in the static case. 
We also obtain the expression for the energy contained within a surface of 
constant radius r Q in the asymptotic region where the metric coefficients may 
be determined. 



Notation: spacetime indices /i, u, ... and local Lorentz indices a, b, ... run from 
to 3. In the 3+1 decomposition latin indices from the middle of the alphabet 
indicate space indices according to \i = Q,i, a = (0), (i). The tetrad field 
e a n and the spin connection u^ab yield the usual definitions of the torsion 
and curvature tensors: R a bpLU = d^uj v a b + c uj u c b - T a ^ = d^e a v + 
^fj, a b£ b v — ■■■■ The flat spacetime metric is fixed by f?(o)(o) — — 1- 
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II. The TEGR 



The Lagrangian density of the TEGR in empty spacetime is given by 



L(e,u,\) = -ke(-T abc T abc + -T abc T bac - T a T a ) + e\ ab ^ R ab ^(u) . (2) 

where k = j^q, G is the gravitational constant; e = det(e a M ), \ ab ^ u are 
Lagrange multipliers and T a is the trace of the torsion tensor defined by 
T a = T b ba . The tetrad field e 0/1 and the spin connection uj^ ab are completely 
independent field variables. The latter is enforced to satisfy the condition of 
zero curvature. Therefore this Lagrangian formulation is in no way similar 
to the usual Palatini formulation, in which the spin connection is related to 
the tetrad field via field equations. Later on we will introduce the tensor T, abc 
defined by 

rpabcrp | rpabcrp rparp rpabc^ 

abc ' 2 1 bac 1 1 a = 1 ^abc ■ 

The equivalence of the TEGR with Einstein's general relativity is based 
on the identity 



eR(e,u) = eR{e) + e(^T abc T abc + l -T abc T bac - T a T a ) - 23 /1 (eT") , (3) 

which is obtained by just substituting the arbitrary spin connection u^ ab = 
°^afe(e) + K^ ab in the scalar curvature tensor R(e,uj) in the left hand side; 
°^Vfc(e) is the Levi-Civita connection and K^ ab = le a x e b u (T XlMU + T vXf i - 
Tfjj/x) is the contorsion tensor. The vanishing of R a b ^ u (u>), which is one of 
the field equations derived from (2), implies the equivalence of the scalar 
curvature R(e), constructed out of e a M only, and the quadratic combination 
of the torsion tensor. It also ensures that the field equation arising from the 
variation of L with respect to e a M is strictly equivalent to Einstein's equations 
in tetrad form. Let = denote the field equations satisfied by e afl . It 
can be shown by explicit calculations that 

= ^{iUe) - \^R(e)} ■ (4) 
We refer the reader to refs. |7|, |8| for additional details. 



4 



Throughout this section we will be interested in asymptoticaly flat space- 
times. The Hamiltonian formulation of the TEGR can be successfully im- 
plemented if we fix the gauge u 0a b = from the outset, since in this case the 
constraints (to be shown below) constitute a first class set@. The condition 
^oab = is achieved by breaking the local Lorentz symmetry of (2). We still 
make use of the residual time dependent gauge symmetry to fix the usual 
time gauge condition en* = emv = 0. Because of u; 0a fc = 0, H does not 
depend on p kab : the momentum canonically conjugated to u> ka b- Therefore 
arbitrary variations of L = pq — H with respect to P kab yields uj ka b = 0. Thus 
in view of u>o a b = 0, u ka b drops out from our considerations. The above gauge 
fixing can be understood as the fixation of a reference frame. 

As a consequence of the above gauge fixing the canonical action integral 
obtained from (2) becomes|| 

A TL = I d A x{U^ k e m - H} , (5) 



H = NC + N l C t + E mn U mn + J—d k (NeT k ) + ^(IF*^-) . (6) 

07TG 

N and N l are the lapse and shift functions, and S mn = — S nm are Lagrange 
multipliers. The constraints are defined by 

C = djpkeTt) - keZ kij T Mj - ^(n^ ~ \^ 2 ) , (7a) 

C k = -e U)k d t U^ 1 - U.WT m , (76) 

with e = det(e(j)k) and T l = g lk e^ l T^i k . We remark that (5) and (6) are 
invariant under global SO (3) and general coordinate transformations. 
If we assume the asymptotic behaviour 

e(i)fe ~ Vjk + \h jk {^) (8) 
for r — > oo, then in view of the relation 

f d 3 xd 3 (eT j ) = f dS k (dih lk - d k h u ) = E ADM (9) 

8ttG J 16nG Js 
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where the surface integral is evaluated for r — > oo, the integral form of the 
Hamiltonian constraint C = may be rewritten as 

J d 3 x{keZ ki iT klJ + -^(n%i - ^n 2 )| = Eadm ■ (io) 

The integration is over the whole three dimensional space. Given that dj(eT^) 
is a scalar density, from (9) and (10) we define the gravitational energy 
density enclosed by a volume V of the space as 



It must be noted that E depends only on the triads etk)i restricted to a three- 
dimensional spacelike hypersurface; the inverse quantities e^ k ' 1 can be written 
in terms of enAi- From the identity (4) we observe that the dynamics of the 
triads does not depend on uj^ ab . Therefore E g given above does not depend 
on the fixation of any gauge for u^. We briefly remark that the reference 
space which defines the zero of energy has been discussed in ref . |§ . 

We make now the important assumption that general form of the canon- 
ical structure of the TEGR is the same for any class of space-times, irrespec- 
tive of the peculiarities of the latter (for the de Sitter spacefill , for example, 
there is an additional term in the Hamiltonian constraint C) . Therefore we 
assert that the integral form of the Hamiltonian constraint equation can be 
written as C = H — E = for any space-time, and that (1) represents the 
gravitational energy for arbitrary space-times with any topology. 

Before closing this section, let us recall that Miiller-Hoissen and Nitschj] 



and Kopczyhski|2T| have shown that in general the theory defined by (2) 
faces difficulties with respect to the Cauchy problem. They have shown that 
in general six components of the torsion tensor are not determined from the 
evolution of the initial data. On the other hand, the constraints of the theory 
constitute a first class set provided we fix the six quantities u 0a b = before 
varying the action^. This condition is mandatory and does not merely 
represent one particular gauge fixing of the theory. Since the fixing of uj 0a b 
yields a well defined theory with first class constraints, we cannot assert that 
the field configurations of the latter are gauge equivalent to configurations 
whose time evolution is not precisely determined. The requirement of local 
SO (3,1) symmetry plus the addition of \ abiiy R a ^ u {uj) in (2) has the ultimate 
effect of discarding the connection. Although we have no proof, we believe 
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that the two properties above (the failure of the Cauchy problem and the 
fixation of u>o a b = 0) are related to each other. 

Constant rotations constitute a basic feature of the teleparallel geometry. 
According to M0ller||, in the framework of the abolute parallelism tetrad 
fields, together with the boundary conditions, uniquely determine a tetrad 
lattice, apart from an arbitrary constant rotation of the tetrads in the lattice. 



III. M0ller's energy expression 

M0ller carried out several investigations regarding the localizability of the 
gravitational energy. He faced difficulties in establishing a covariant expres- 
sion using the metric tensor 0, and because of this he arrived at an expression 
through the use of tetrads @, [|. According to M0ller, this latter expression 
still has a difficulty in that it is not invariant under local Lorentz trans- 
formations. It is very instructive to compare expression (1) with M0ller's 
expression. For the sake of this comparison, we will put aside the difficulty 
regarding the noninvariance with respect to local Lorentz transformations. 

M0ller presents an expression for the energy-momentum of the gravita- 
tional field. However we will only consider the energy expression. Translating 
into our notation, Miller's energy reads 

E = -Jd 3 xd x U 0X , (11) 
where the potential in the integrand is given by 

U, vX = [e™V»e a A + (Sy x - 8^) V CT e a CT ] (12) . 

In contrast with the notation of the previous section, all geometrical quanti- 
ties in equations (11-15) are four- dimensional quantities. In (12) V repre- 
sents the covariant derivative with respect to the Christoffel symbols F* v . 
M0ller's energy can be first rewritten as 

E = f d^xd^ee^V^) . (13) 

By means of the identity 
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V fc e aj = d k e aj - Tljeaa = -°u k a e bj , (14) 
where "uj^ab is the Levi-Civita connection, 

^fjiab 2^ /^(^afcc ^bac ^cab) ; 

^abc ^-avifib^ d^6 c G-c^d^Gfr ) , 

we can further rewrite expression (13) as 

E = -^qJ ^xdiiee^^jab) . (15) 

Up to this point {e aAt } are tetrads of the four- dimensional space-time. In 
order to compare (15) with (1) let us impose the time gauge f, = e^) = 
and establish the 3+1 decomposition of the tetrads as in ]7|, 22]. Then the 
integrand on the right hand side of (15) can be rewritten as 



4 e A e ai A e bj ° 



u jab ( 4 e) = Nee^e™ °u j{m){n) (e) - e{e ai W - e a > N l ) °u m{m) 



N and N l are the lapse and shift functions and { 4 e af1 } are tetrads of the four- 
dimensional space-time. In terms of triads restricted to a three-dimensional 
spacelike surface we have 

E = 8^G I d3xd ^ Nee{m)lein)3 °^M(n) ~ e(e (m)i i^ - e ^JV*) o w i(0)M ] . 

(16a) 

Comparison with (1) can now be made if we make use of the identity 



d^ee^e^ °u j{m){n) ) = d t (eT) , 

where the right hand side above is the same as in (1). M0ller energy can be 
finally written as 

E = 8^G I d * xd ^ NeTi ~ <z [m)lN3 ~ e {m)j N l ) °cu m{m) )} ■ (166) 
Recall that we are ignoring local Lorentz transformations. 
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Besides the appearance of extra terms involving °Wj(o)(m) 011 the right 
hand side of (16a,b), there is also the crucial presence of the lapse function 
N multiplying eT\ Therefore even for configurations of the gravitational field 
for which the second term on the right hand side of (16) does not contribute 
(if, say, N l = 0, as for the Schwarzschild solution) expressions (1) and (16) 
will lead to different results when applied to finite volumes of the three- 
dimensional space. Moreover, because of the presence of the lapse function, 
(16) is not invariant under time reparametrizations: N'(x'°) = ^q-N(x°). 
Thus for a finite volume of integration (16b) does not remain invariant under 
this reparametrization. 

In the Einstein-Cartan theory the connection °u;j(o)(m) can be expressed 
in terms of the momenta canonically conjugated to er m )i- in the notation of 
2|] it is given by °a>j(o)( m ) = h( n (m)j ~ \^{m)j^) (see equation 12 of p2| ). In 



this context (16b) can be rewritten as 



E = d^NeV - ieW% Wj ] . 

The expression above is exactly the energy expression for the Einstein-Cartan 
theory (see equation 21 of f22|), assuming that the gravitational energy is 



obtained from integration of surface terms of the total Hamiltonian. This 
expression is also very similar to (i) the integral of the surface terms in 
equation 6, and (ii) the energy expression considered by Nester[^|] in the 
analysis of the positivity of the gravitational energy (equation (3.15) of p4| ). 
All definitions of gravitational energy considered above agree for the total 
gravitational energy. 



IV. Bondi's radiating metric and the associated triads 

Bondi's metric is a not an exact solution of Einstein's equations. It de- 
scribes the asymptotic form of a radiating solution. In terms of radiation 
coordinates (u, r, 9, <fi), where u is the retarded time and r is the luminosity 
distance, Bondi's radiating metric is written as 

ds 2 = -(^-e 2f3 -U 2 r 2 e 2 ^du 2 -2e 2f3 dudr-2Ur 2 e 2 ^dud9 
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+r 2 (e 2 ^ d9 2 + e~ 2 ~< sin 2 9 d^J . (17) 

The metric above is such that the surfaces for which u = constant are null 
hypersurfaces. Each null radial (light) ray is labelled by particular values of 
u, 9 and <p. At spacelike infinity u takes the standard form u = t — r. The 
four quantities appearing in (17), V, U,(3 and 7 are functions of u, r and 9. 
Thus (17) displays axial symmetry. A more general form of this metric has 
been given by Sachs p5fl , who showed that the most general metric tensor 
describing asymptotically flat gravitational waves depends on six functions 
of the coordinates. 

The functions in (17) satisfy the following asymptotic behaviour: 



P= -S3 + - 



7 = ; + 0(^) + ... 



7 = 1 



2M 1 



dd /<9c\ 2 fdc\ If 

_ +d coW - y -ic( w )cow-- <> (i+scoee 



+ . 



1 /(9c 
C/ = — ( — 4 2c cold 



1 / Be 
+ — f2d + 3c — + 4c 2 col 9 



+ 



.00 / r 3 V 50 

where M = M(tt, 9) and ci = d(-u, 9) are the mass aspect and the dipole 
aspect, respectively, and from the function c(u, 9) we define the news function 

dc(u,B) 
du 

The application of (1) to Bondi's metric requires transforming it to co- 
ordinates t, r, 9 and 4> for which t = constant defines a space-like hypersur- 
face. Before proceeding, we recall that the analysis of (17) in t,x,y,z co- 
ordinates has already been performed by Goldberg[^], in the investigation 
of the asymptotic invariants of gravitational radiation fields. Therefore we 
carry out a coordinate transformation such that the new timelike coordinate 
is given by t = u + r. We arrive at 



ds 2 = -(Te 2 ^-U 2 r 2 e 2 ^yt 2 -2Ur 2 e 2 ^dtd9+2^(^-l^ 



-U 2 r 2 e 2 ^ 



dr dt 
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+ 



V 



e 2/3 2--)+U 2 r 2 e^ 
\ r 



dr 2 + 2U r 2 e^dr d9 + r 2 [e^dO 2 + e" 27 sin 2 6 dcf) 2 ) . 

(18) 

Therefore the metric restricted to a three-dimensional spacelike hypersurface 
is given by 



ds 2 



o2/3 



2-^) +U 2 r 2 e 2 ^ 
r 



dr 2 + 2Ur 2 e 27 dr d9 



(19) 



We must consider triads that correspond to the metric above. The con- 
struction of triads, in general, is a nontrivial step. If in a given coordinate 
system the metric tensor is diagonal, then the construction of triads is a rel- 
atively simple procedure. One must only make sure that the triads satisfy 
the appropriate boundary conditions at infinity. Recall that in order to ob- 
tain expression (9) for the ADM energy the triads must have the appropriate 
asymptotic behaviour given by equation (8). 

The metric tensor (19) has an off-diagonal element, and this fact adds 
a bit of complication in the construction of triads. Nevertheless we can 
immediately write down two sets of triads that lead to this metric. They are 
given by 



e(k)i 



where 



( A sinO coscf) + B cos8 coscf) rC cosO coscj) —rD sinO sin(f> \ 
A sinO sincj) + B cosO sine/) rC cosd sincj) rD sinO coscf) 
A cosO — B sinO —rC sinO 



(20) 



A = e\h- 



B = rUe^ 



V 



(21a) 
(216) 



C = e 7 , 



(21c) 
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D 



(21c/) 



and 



where 



/ A' sin9 cose, 
A' sinO sine 
\ A' cos9 



rB' cos9 coscj) + rC sin9 cose, 
rB' cos9 sincj) + rC sin9 sine, 
—rB' sin9 + rC cos9 



—rD' sin9 sine] 
rD' sin9 coscj) 




A' 



+ U 2 r 2 ^ 




A' 



D' 



(22) 

(23a) 

(236) 

(23c) 
(23d) 



It is easy to see that both (20) and (22) yield the metric tensor (19) through 
the relation e^e^fc = g^. They are related by a local SO (3) transformation. 

Triads given by (20) and (22) are the simplest sets of triads that satisfy the 
two basic requirements: (i) the triads must have the asymptotic behaviour 
given by (8); (ii) by making the physical parameters of the metric vanish we 



must have Ti 



(k)ij 



everywhere. In the present case if we make M = d 



both (20) and (22) acquire the form 



e(fc)i 



/ sin9 cose/ 
sin9 sine, 

V 



cost 



r cos9 cose, 
r cos9 sine, 
—r sin9 



— r sin9 sine, 
r sin9 coscj) 




(24) 



In cartesian coordinates the expression above can be reduced to the diagonal 
form e(fc)j(x, y, z) = The requirement (ii) above is essentialy equivalent 
to the establishment of a reference space, as discussed in ||. Without the 
notion of a reference space we cannot define gravitational energy. Note that 
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by a suitable choice of a local SO (3) rotation we can make the flat space 
triad (24) satisfy the requirement (i), but not (ii). 

It is possible to show that (20) and (22) are particular cases of an infinit 
set of triads that satisfy both requirements above. It is given by 



/ A sinOcoscf) + B cos9coscf> rC cos9coscf> + rT> sinOcoscf) — re~ 7 sinOsincJ) 
AsinOsincj) + B cosOsincf) rC cosOsintp + rT> sinOsincf) re -1 sin9cos(j) 
\ A cosO — B sinO —rC sinO + rV> cos9 



e(fc)t 

with the following definitions: 



(25) 



A = J&^-'jj +U 2 r 2 (e 2 'f -b 2 
B = bUr , 
C = Ve 2 ^ - d 2 U 2 r 2 , 
V = dUr , 

where b and d are arbitrary, dimensionless functions that must satisfy 



V 



b = Ve^-d 2 U 2 r 2 + de" 7+/3 y2 

By making d = we obtain (20), and b = leads to (22). 

From the point of view of the TEGR triads given by (20) and (22) are 
physically inequivalent (that is, they are not gauge equivalent), because we 
have seen that the Hamiltonian formulation established by eq.(5) is not in- 
variant under the local SO (3) group, but rather under the global SO (3). In 
the TEGR the torsion tensor describes the way in which the space-time is de- 
formed. The latter is thus considered as a continuum with microstructure[5|]. 
Therefore the same space, defined uniquely by the metric tensor, may be de- 
formed in several ways, according to the manner one defines the triads. This 
is essentially the geometrical meaning of the noninvariance of the TEGR 
under local SO (3) transformations. 



13 



In the Hamiltonian formulation of the TEGR the basic geometrical field 
variable is the triad, not the metric tensor. Any set of triads should be ruled 
out on physical grounds, i.e., if they lead to incorrect physical statements 
concerning the energy content of the gravitational field. 

In the next section we will obtain the expressions for the gravitational 
energy arising from (20) and (22). These expressions are quite different. 
Although the expression corresponding to (20) is simpler, as we will see, we 
have no definite experimental evidence in favour of it. 

V. Gravitational radiation energy 

In this section we will apply expression (1) both to (20) and (22). Since 
the two triads display distinct geometrical properties, we expect to obtain 
different expressions for the energy density -^di(eT l ) (we will make the grav- 
itational constant G — 1). Our analysis is meaningful only in the asymptotic 
region of large values of the radial distance. However, we have no reason to 
expect (20) and (22) to yield different expressions for the total energy of the 
field. In fact, as we will see, they yield the same (expected) expression. 

As we mentioned earlier, the significance of the present approach to the 
analysis of gravitational radiation fields resides in the fact that we can eval- 
uate the gravitational energy inside a large but finite portion of a three- 
dimensional spacelike surface. In other words, by means of Gauss law (1) can 
be evaluated over a surface far from the source of radiation, in the asymptotic 
limit where the metric components are precisely determined. Specifically, we 
will evaluate (1) inside a large sphere of radius r Q . The time evolution of the 
metric field will determine the time dependence of this energy, and conse- 
quently the energy radiated out of it. 

In order to calculate T % = g lk 'g m: > 'e^ m T^j k we need the inverse metric g % K 
In terms of the definitions (21), it is given by 



We will initially consider the set of triads given by (20). In the following, 
a comma after a field quantity indicates a derivative: A,i and A, 2 indicate 




(26) 
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derivative with of A respect to r and 9, respectively. The torsion components 
for (20) are given by 



^(1)12 = cos 9 coscj) (C + rC,i —A — B, 2 ) + sinO coscj) (— A, 2 +B) 



^(1)13 — sin® sincj) (A — D — rD,i ) + cos9 sincj) B , 
2~(i)23 = sin9 sincf) (— rD, 2 ) + cos9 sincj) (rC — rD) , 



T(2)i2 = cos9 sincp (C + rC,\ —A — B, 2 ) + sin9 sincf) (—A, 2 +B) , 



T( 2 )i3 = sin9 coscj) {—A + D + rD,i ) + cos9 coscj) (-B) , 
T(2)23 = sin9 coscj) (rD, 2 ) + cos9 coscj> (— rC + rD) , 



7(3)12 = sin0 (A + B, 2 -C - rC* ) + cos9 (-A, 2 +B) , 



T(3)13 = ^(3)23 = . 

Since we are interested in calculating the energy inside a surface of con- 
stant radius, only T 1 will be considered. By Gauss law, the expression of 
this energy is given by 



E = — [ d9dcj>eT 1 , 

87T Js 



(27) 



where S is a surface of fixed radius r D , assumed to be large as compared with 
the dimension of the source, and the determinant e is given by e = r 2 A sin9. 
After a long but otherwise straightforward calculation we arrive at 



Ej = -± / d9 I sm9 

4 jo 



e 7 + e 



_ 7 _ 



A 



+ j^(Ursin9)\ . (28) 



with A defined by (21a). Let us note that the field quantities appearing in 
(28) are functions of u — t—r: M = M(t—r, 9),c — c(t—r, 9), etc. Therefore, 
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once these functions are known, one can explicitely calculate the variation 
of Ei with respect to the time t, however only in the limit where the metric 
components are precisely determined. 

Unfortunately the expansion of Ei up to terms in — , making use of the 
aymptotic behaviour of U, V, (3 and 7, yields no simple expression. It is given 
by 



Ej 



- r fflsinOM - — f ' d6 sin9 
2 Jo 4r c Jo 



1 

4r, 



/"7T 

o JO 



d9M 



d_ 

89 



sinO 



8c 
89 



2ccot9 



(29) 



In the calculation above we have assumed that U(9)sin9 = d{9)sin9 = 
when 9 = or n. 

We observe that Ei yields Bondi energy in the limit r — > 00 in the static 
case (i.e., when M is a function of 9 only), so that E t is the total gravitational 
energy. However, in the nonstatic case an expression for the loss of mass due 
to gravitational radiation can be obtained from (29). This is one major result 
of our analysis. 

Let us recall that Bondi's mass aspect m(u), 



m(u) = \ I d9sin9 M(u,9) , 
2 Jo 



(30) 



depends on the null foliation used. The mass aspect m(u) can be understood 
as a mass associated to each null cone determined by the equation u = 
constant. Since the limit r — > 00 corresponds to«^ —00 for finite t, we see 
once again that in this limit Ei gives the total energy because it corresponds 
to the initial value of the Bondi energy. 

We will consider next the second set of triads, eq.(22). The components 
of the torsion tensor resulting from the latter are given by 



T (1)12 = cos9 cos<p (B 1 + rB',i -A') + sin9 cos<p (C + rC ',1 -A', 2 ) , 



7(1)13 = sinO sin(i> {A' — D' — rD' :1 ) , 
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2^(1)23 — cos & sincp (rB' — rD') + sinO sine/) (rC — rD' , 2 ) , 



T (2)12 = cosd sine/) (B' + rB' a -A') + sinO sirnf) (C + rC",i -A', 2 ) , 

T (2)13 = -sinO coscp {A' — D' — rD' a ) , 
T (2)23 = -cos9 cos</> (rB' - rD') - sine coscp (rC - rD', 2 ) , 
7(3)i2 = sine (B' + rB' a -A') + cose (C + rC" a -A', 2 ) , 

7(3)13 = 7(3)23 = . 

As in the previous case, we are interested in calculating the energy in the 
interior of a surface of constant radius r Q . Therefore only the knowledge of 
T 1 will be necessary. After a long calculation we first arrive at 



eT 1 = 



r sint 



A 



-27 



+ e 



rB sine 



27 



c?7 
dr 



r (9-7 


rB cose 






B' - e" 7 




A 





(31) 



where the primed quantities are given by (23). It is not straightforward to 
put the expression above in a simplified form. After some rearrangements 
we can finally write the energy expression (27) as 



E 



11 



sine 



-27 m 



r A' + e-^A'B' -2 + e-^—C - BC' - Bc^< — '- 



,<97 



69 



-B cose 



B' - e" 7 



de 

(32) 
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Like equation (28), this expression represents the energy enclosed by a large 
spherical surface of radius r Q . Expanding the expression above up to the first 
power of — we find 



E 



11 



1 r 

2 Jo 



dd M sin6 / ddsinO 

4r G Jo 



3M 2 + 



5 fdc\ 
2\d9~) 



1 2 • L0c(|f]< 



89 



+8c 2 cot 2 e - ( ( ^ + 2c cote] 1 - — r de cose + 4c 2 cote 

V 09 J \a9 J J Ar Q Jo i \o6 J 

(33) 

We are again assuming U(e)sin,e = d(8)sine = for 9 = 0, ir. 

We observe that in the limit r — > oo E 1 // also gives the total energy 
As before, for a finite (but sufficiently large) value of r G we can compute the 
loss of mass due to gravitational radiation, once the functions M and c are 
known in the asymptotic region. 



VI. The selection of triads 



In section IV we obtained an infinit set of triads that yield the three- 
dimensional spacelike section of Bondi's metric, and in the previous section 
we considered in detail only the simplest constructions. Of course simplicity 
is a major feature of physical systems, but we are really in need of experimen- 
tal evidence that leads to a definite description. We need actual realizations 
of the quantities M(r — t,6), c(r — t,9), d(r — t, 9) and experimental evidence 
on how the energy is radiated away in order to arrive simultaneously at the 
correct energy expression arising from (27) and at the definite expression of 
e(k)i- 

However, we can envisage two possible types of conditions on the triads 
that associate a unique triad with the three-dimensional metric tensor. 

The first condition regards the energy content of the gravitational field. 
If we stick to the point of view according to which physical systems in nature 
have a tendency to be in states of minimum energy, then the correct triad for 
the spacelike section of Bondi's metric is the one that minimizes expression 
(27) for all possible constructions of e^)i- By means of this criterium we 



18 



consider triads given by (25), or any further construction that complies with 
the two conditions stated in section IV, and ask which one yields the smaller 
value of energy contained within a surface of constant radius, in similarity 
with the calculations of (29) and (33). Unfortunately this analysis cannot be 
carried out unless M, c and d are known. 

Certainly one can ask whether only (27) should be minimized or the 
energy density should be everywhere a minimum. In the context of Bondi's 
metric the latter possibility cannot be considered, because the metric is valid 
only in the asymptotic limit, but in the general case it is an open question 
that must be carefully addressed. 

The second condition takes into account equation (8): we require the 
triads to have the asymptotic behaviour determined by (8) with a symmetric 
tensor hj k = h k j. Again, one has to find out of (25) which realization of 
e(j)fe in cartesian coordinates complies with this criterium. This condition 
may be understood as a rotational gauge condition. Note that as it stands, 
hjk in equation (8) is not required to be symmetric (in equation (9) only 
the symmetrical part contributes). By explicit calculations we observe that 
neither (20) nor (22) satisfy this second condition. 

The two conditions above may not be mutually excluding. On the con- 
trary, they may lead to the same triad. The determination of the correct 
triad is certainly an essential and crucial issue of the theory and will be 
further investigated elsewhere in the general case, with special attention to 
Bondi's metric, in the light of the conditions above. 

We observed that both (29) and (33) yield the same total energy. This 
is also the case if we carry out the calculations with a more complicated 
triad, whether belonging to (25) or not, which is related to (20) or (22) by a 
local SO (3) transformation with an appropriate asymptotic behaviour. Let 
us consider a local SO (3) transformation given by 



Expression (35) can be best analysed if we consider an infinitesimal rotation. 




(34) 



E = E + — J^xdAeg lk ^ l \ m) e {m)] {e [n)k d 3 K {l) 



(n) 



e (n)j d k A {l) (»>)] . (35) 
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We assume that in the limit r — > oo the SO (3) elements have the asymptotic 
behaviour 



such that 0,1 U(k)(i) = — °' 1 LJ(i)(k), and { cu^)(i)} are constants. Taking into 
account (8) it is easy to see that when integrated over the whole three- 
dimensional space the integral on the right hand side of (35) reduces to a 
vanishing expression: 



Therefore we expect to find the same result for the total energy if we evaluate 
(27) out of any element of the set of triads (25). 

VII. Discussion 

The application of the energy definition (1) for a given solution of Ein- 
stein's equations requires considering a foliation of the space-time in three- 
dimensional spacelike surfaces. The metric for the spacelike section of Bondi's 
radiating metric admits an infinit set of triads related by local SO (3) trans- 
formations. In the present case, from this infinit set of triads we singled out 
two of them. We have considered in detail the two ones that exhibit the 
simplest structures in spherical coordinates, and that: (i) satisfy the asymp- 
totic conditions given by (8); (ii) reduce to the reference space {T(k)ij — 
everywhere) if we make the physical parameters vanish: M = c = d = 0. 

The two sets of triads, (20) and (22), describe the spacelike section of 
Bondi's metric given by (19) and lead to the energy expressions (28) and 
(32), respectively. These expressions establish distinct and quite definite 
physical predictions. They allow us to compute the energy radiated from the 
interior of a spherical surface of constant radius r . It would be a remarkable 
achievement of the TEGR if, on physical grounds, we could decide for one of 
them or even for an arbitrary element of (25). In the TEGR the space (space- 
time) geometry is fundamentally described by triads (tetrads). Unfortunately 
we do not dispose of experimental information for taking such a decision. 

It is a very important result that the total energy due to both sets of triads 
(as well as from any element of (25) ) agrees exactly with the static Bondi 
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energy, in which case the energy arises from the integration of M — M{9). 
In fact, the definition of Bondi's mass aspect is basically motivated by the 
fact that in the static case, and by investigating the asymptotic properties 
of the gravitational field, M{6) arises as the mass of an isolated system. 

The final expressions (29) and (33) support the consistency of the defi- 
nition (1), and the relevance of the TEGR as a fundamental description of 
general relativity. However the present analysis, which was developed on 
spacelike surfaces, has to be compared with the one recently carried out on 
null surfaces |27|. Let us recall that in order to obtain the Hamiltonian for- 
mulation given by (5-7) we imposed the time gauge condition. Therefore the 
resulting geometry may be understood as a three-dimensional teleparallel ge- 
ometry, since the teleparallelism is restricted to the three-dimensional space- 
like surface. On the other hand, in the Hamiltonian formulation developed 
in [E7| we have not fixed any particular tetrad component, and consequently 
the teleparallel geometry is really four- dimensional. 

In p7j the constraints also contain a total divergence, in similarity with 
(7a), and may be taken likewise to define the gravitational energy-momentum. 
Although it appears that the geometrical framework of |27| is better suited to 
the analysis of the Bondi-Sachs metric, we note that the the energy expres- 
sion arising there is considerably more complicated than (1). Moreover we do 
not know yet whether the constraint algebra leads to a consistent Hamilto- 
nian formulation (either on null or spacelike surfaces). We also note that one 
Hamiltonian formulation cannot be reduced to the other by means of gauge 
fixing. Nature admits only one correct physical description, and therefore 
either the three-dimensional or the four- dimensional teleparallel geometry is 
the correct candidate for describing the energy properties of the gravitational 
field. All these issues will be considered in the near future. 
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